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An  Analytic  Theory  for  Trajectories  and  Current 
to  a  Cylinder  in  a  Flowing  Magnetoplasma 


I.  INTRODUCTION 


Accurate  comprehensive  equations  for  the  electrical  behavior  of  spacecraft 
and  plasma  in  the  vicinity  cannot  be  solved  analytically.  Therefore,  either  they 
have  been  solved  numerically  with  time  consuming  computer  programs,  or  many 
simplifying  assumptions  have  been  made. 

The  theory  of  electric  probes  with  simple  geometry  in  a  stationary  plasma 
with  or  without  a  magnetic  field  has  been  treated  fairly  fully.  2,  3,  4,  5,  6 
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Less  work  has 


J  .  _  .  ,  7.8.9,10.11 

work  has  also  been  done  on  probes  in  flowing  plasmas. 

been  done  on  the  more  complicated  problem  of  a  probe  in  a  flowing  magnetoplasma. 

Rough  calculations  of  the  effect  of  magnetic  field  B  on  satellite  drag  were  made  by 

Brundin^  and  Drell  et  al.  ^  Whipple  ^  showed  that,  as  the  gyroradius  decreases 

from  100  to  0.01  of  the  sphere  radius,  the  current  to  a  sphere  decreases  to  one- 

half.  By  numerical  solution  of  Vlasov's  and  Poisson's  equations,  Grabowski  and 
15  16 

Fischer  °  and  Brooks  and  Koehler  °  obtained  the  ion  density  dependence  on  radius 
and  azimuth  of  a  cylinder  parallel  to  the  magnetic  field  moving  perpendicular  to  it. 
Grabowski  and  Fischer  found  that,  for  equal  spacecraft  and  ion  thermal  speeds, 
the  magnetic  field  became  quite  important  when  the  gyroradius  decreased  to  the 
cylinder  radius.  They  did  not  consider  the  case  of  smaller  gyroradii.  Brooks  and 
Koehler  considered  sounding  rockets  with  radius  smaller  than  the  ion  gyroradius 
and  so  neglected  magnetic  field  effects. 

In  this  paper,  simplifying  assumptions  are  made  and  a  new  method  explored 
that  permits  an  analytical  solution  of  the  electrical  characteristics.  Application  to 
low  equatorial  Earth  satellites  larger  than  the  shuttle  is  envisioned.  The  basic 
theory  is  presented  in  the  next  section.  It  is  then  applied  to  a  cylinder  parallel  to 
the  magnetic  field. 
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2.  BASIC  THEORY 


This  paper  concerns  the  case  of  a  conducting  body  moving  perpendicular  to  a 
known  uniform  magnetic  field  with  constant  velocity  VQ  in  a  magnetoplasma.  The 
formalism  for  describing  the  motion  of  the  plasma  particles  is  the  guiding  center 
approximation.  This  is  valid  when  the  gyroradius  times  each  component  of  elec¬ 
tric  field  gradient  is  much  smaller  than  the  electric  field.  Except  where  the  elec¬ 
tric  field  nearly  vanishes,  this  holds  when  the  gyroradius  is  small  compared  to 
the  body  dimensions.  The  effect  of  the  finite  size  of  the  gyroradius  is  neglected. 
As  in  Northrup,  the  lowest  ("zeroth")  and  next  ("first,"  m/e  times  zeroth)  order 
terms  are  kept,  and  particle  motion  is  averaged  over  a  gyroperiod.  Steady  lami¬ 
nar  flow  is  assumed,  and  plasma  oscillations  are  neglected.  Boltzmann's  equation 
is  not  used.  The  ion  thermal  speed  is  smaller  than  the  ion  drift  speed  and,  for 
simplicity,  is  neglected.  The  treatment  is  two-dimensional,  and  thus,  for  exam¬ 
ple,  is  applicable  to  the  magnetic  field  parallel  to  a  cylinder  of  length  large  com¬ 
pared  to  other  dimensions.  The  ambient  plasma  is  assumed  to  be  uniform,  colli¬ 
sionless,  and  without  electric  field.  The  magnetic  and  the  electric  fields  are  as¬ 
sumed  to  be  time  independent  in  the  frame  of  the  body,  the  frame  in  which  the  cal¬ 
culation  is  made.  Therefore,  the  electric  field,  E  =  -  V4>. 

Only  nonrelativistic  velocities  are  considered;  the  magnetic  field  caused  by 
currents  is  negligible  and  so  is  neglected.  All  charges  reaching  the  body  are  as¬ 
sumed  to  be  absorbed;  no  scattering,  photo,  or  secondary  emission  is  assumed. 
The  latter  assumption  is  valid  for  the  low  particle  energies  considered  here.  For 
simplicity,  the  force  of  gravity  is  neglected  since  it  is  relatively  small.  In  MKS 

units  with  the  z-axis  chosen  in  the  direction  of  the  magnetic  field  B,  the  drift  ve- 
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locity  of  the  center  of  gyration  of  a  particle  is  given  by  Eq.  (1.  17)  of  Northrop 


Vi 


M 

e 


m  dz 

FB  *  —  (-g  +  v|(—  + 
e  dt 


(1) 


M,  m,  e,  and  V|j  are  the  particle's  magnetic  moment,  mass,  charge  (including 
sign),  and  component  of  velocity  parallel  to  B  respectively;  ~g  is  the  acceleration 
of  gravity. 


17.  Northrop,  T.G.  (1963)  The  Adiabatic  Motion  of  Charged  Particles, 
Interscience  Publishers,  p.  S.  ~  * 


i 


v  -  E  *  B 
E~  B2 


is  the  zeroth  order  drift  velocity.  From  the  assumptions  above,  the  second,  third, 
and  fourth  (middle)  terms  of  the  right-hand  side  of  Eq.  (1)  are  neglected  and  V  » 
V^.  Since  B  is  constant  and  d/di  -  0,  the  last  term  becomes  m/(eB^)V  •  7  E. 

So  (to  the  first  order),  Eq.  (1)  becomes 


V  i  v  +  V 

E  1* 


where 


v"  =  — rn  "v  •  7  eT 
1  eB2  E 


The  latter  is  the  first  order  drift  velocity,  often  called  the  "polarization"  drift. 

At  the  surface  of  the  conductor,  E  is  perpendicular  to  the  conductor,  so  Vj_,  is  tan¬ 
gential  to  the  conductor.  Therefore,  to  zeroth  order,  the  charged  particles  go 
around  the  conductor  without  hitting  it.  The  total  guiding  center  energy,  £ ,  of  a 
particle  is  a  constant  of  motion.  Since,  to  first  order. 


£  -  — v  +  e<*. 
2  E  v 


this  is  constant  along  the  guiding  center  trajectory,  that  is,  along  a  streamline  of 
the  specie.  For  this  two-dimensional  case. 


V  =  _i_  B  x  7e 

i  m  r7 

eH“  1  - - v  •  y 

eH" 


is  consistent  with  F'q.  (3).  This  may  be  seen  by  writing  out  the  four  terms  in 

2 

Eq.  (4)  from  the  components  of  the  vectors,  substituting  1  «•  m/(eR  )  V  •  E  for 
[l  -  m/(eB“)7  •  Ej  in  Eq.  (6),  using  Eq.  (2)  in  Eq.  (5),  which  is  then  substi¬ 
tuted  into  Eq.  (6),  and  writing  out  the  terms.  From  Eqs.  (3),  (4),  and  (6),  the 
current  density  of  specie  or  is 


na  ca  va 


nnca  ^  •  +  <Vr,*  7)E 


'  NY.'/*’ 


mm  mm  warn  mm  msm  af  msm  isaagsa  mm 


where  na  is  the  number  density 


Na  = 


1  - 


ma 


V»E 


ea® 


(9) 


and  the  mass  density  is 
Pa  =  mana- 


(10) 


Taking  the  divergence  of  Eq.  (8)  and  using  vector  identities, 

ea^  *  ^nava  ^  =  — o  *  ^  ^ 

B2 


(11) 


From  charge  conservation  (the  continuity  equation),  Eq.  (11)  =  0.  From  this,  Na 
can  be  proven  to  equal  Na(ea),  that  is,  to  depend  on  coordinates  only  through  its 
dependence  on  Ea.  Thus,  it  is  a  constant  along  a  streamline.  Since  Na  =  na  at  <» 
it  is  independent  of  ta  and  is  the  same  (a  constant)  everywhere  that  the  streamline 
goes,  that  is,  everywhere  except  in  forbidden  regions  and  the  part  of  the  wake 
that  is  free  from  these  streamlines  because  of  sweptout  impacting  particles. 

2 

Writing  Eq.  (9)  as  Na  =  na[l  '■(ma/eaBi^/eo]  from  expansion  and  Poisson's 
equation,  multiplying  it  by  ea  and  summing  over  a. 


-  7^5’- 


(12) 


where  t?  =  and  p  =  £  p  .  Since  (1)  the  left  side  of  Eq.  (12)  is  the  same 

a  u  u  a  a 

everywhere  except  in  the  regions  in  which  streamlines  of  one  or  more  species  are 
missing,  <2l  it  vanishes  at  »  and  thus  also  in  this  region,  and  (3)  the  parenthesis 
on  the  right-hand  side  is  always  positive,  then  rj  must  vanish  throughout  this  re¬ 
gion.  Therefore,  V  •  F.  0  from  Poisson's  equation,  and  Na  -  na  from  F.q.  (9). 
The  current  to  Hie  surface  of  the  bodv  is 


where 


from  Eqs.  (2)  and  (7).  A  is  the  surface  area  of  the  body,  and  E  is  the  electric 
field  at  the  surface.  Jan  is  seen  to  be  independent  of  ea  and  to  come  only  from 
the  polarization  drift  (first  order,  no  2eroth  order  part).  Thus,  the  charging  cur¬ 
rent,  I,  also  is  only  a  first  order  current.  Let  ~y  be  the  distance  vector  from  the 
center  of  the  spacecraft  in  the  B  x  VQ  direction.  If,  for  each  specie  a,  all  and 
only  particles  upstream  having  yam  <  yQ  <  at  —  bit  the  body,  then  it  may 

be  simpler  to  obtain  the  current  per  unit  length  from 


!L  =  a§Ja«<yaM  "  ^am’’ 


where  Jo-  is  the  current  density  far  upstream. 


(15) 


•t.  XIM’I.IC  ATIOM 

The  theory  is  applied  to  a  right,  circular,  conducting  cylinder  of  radius  a  and 
half  length  ba.  (All  lengths  except  a  and  Rg  are  normalized  to  a. )  The  origin  of 
the  coordinate  system  is  taken  at  the  center  of  the  cylinder,  the  x-axis  in  the  di¬ 
rection  of  the  flow  velocity,  -VQ,  in  the  cylinder  frame,  and  the  z-axis  in  the  B 
and  cylinder  axis  direction  (Figure  1).  First,  the  potential,  ,  then  the  energy 
per  particle,  E  ,  are  calculated  in  the  cylinder  frame.  The  potential  resulting 
from  nonzero  charge  density  in  the  sheath,  the  forbidden  region  (see  below),  and 
in  the  part  of  the  wake  with  missing  streamlines  resulting  from  ram  sweepout  is 
neglected,  although  it  must  be  significant.  The  potential  is 

<f>  =  0  4  0,  (16) 

p  q 

where  0p  is  that  due  to  the  electric  field  from  the  motion  and  the  polarization 

charge  on  the  cylinder  caused  by  this  field,  and  0^  is  the  potential  caused  by  a  net 

charge  q  on  the  cylinder.  $  is  the  solution  of  Laplace's  equation  with  z  =  0  for 

P  ^ 

the  cylinder  with  no  net  charge  in  a  background  electric  field  VQBy: 


</>  -V  BRaf  sin  <t>. 


(17) 


PLASMA 


/  PLASMA 


B  B 

a.  stationary  plasma  b.  stationary  cylinder 

Figure  1.  Relative  Motion  Between  Cylinder  and  Magnetoplasma:  Two  Frames  of 
Reference 


where  <t>  =  ^  x,  R.  For  R<<  b,  f  =  1  -  1  /  R 2 ;  for  R  >>  b,  f  ■  1  -  b/R3.  For 

p 

all  ranges  of  R.  f  is  taken  as  p  -  b/[R  (b  +  R)],  where  p  =  b/(b  +  1)  (to  make  A 
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vanish  for  R  =  1).  An  empirical  formula  for  z  =  0  is  used  for 

a  =c-IS[<i  +  |)u-i].  da: 

h  a  n 

where  g  =  1/  (4rr<o),  c  =  1.  777,  and  u  =  l/(cb).  So  the  potential  is  taken  as 


18.  Dubs,  C.  W.  (1982)  Potentials  and  Charges  on  Conducting  Rocket  Sections, 
AFG L-TR -82 -0349,  AD  A 130 143,  p.  29,  Eq.  (All). 


Figure  2.  Equipotential  Surfaces  for  #c  =  0.  <P  -  0  on  the 
plane  with  the  missing  contour 


Figure  3  is  the  same,  except  that  0c  >  0.  0  <  0  for  y  greater  than  the  value 
on  the  missing  contour,  and  0  >  0  for  y  smaller.  Because  of  the  small  mass  of 
electrons,  Vj  <  <  VE>  so  ambient  electrons  nearly  follow  the  equipotential 
curves. 


Y 


X 


Figure  3.  Equipotential  Surfaces  for  0C  >  0.  0  -  0  on 

the  open  cylindrical  surface  with  the  missing  contour 


From  Eqs.  (2)  and  (5), 
£  m(E2  +  E^) 

~  =  - 5—^  +  <P- 

e  2eBz 


(22) 


Since  £  is  the  same  if  0  be  replaced  with  tt  -  0,  it  and  the  trajectories  that  miss 
the  cylinder  are  symmetric  in  x.  This  means  that  the  charge  distribution  in  the 
wake  is  the  mirror  image  of  that  in  the  ram  (the  x  and  y  components  of  velocity  at 


f: 


symmetric  points  are  the  same  except  for  the  sign  of  the  y  component)  except  that 
the  wake  is  missing  the  particles  that  hit  the  cylinder.  Ion  focusing  in  the  wake  is 
not  obtained  because  of  omitting  the  potential  of  the  electrons  in  the  forbidden  re¬ 
gion  and  wake.  For  VD  *  8km/sec,  y  ■  0,  and  R  =  •»,  E/e  is  0.00018V  for  elec¬ 
trons  and  5.3  V  for  0+  ions.  The  gyroradius  of  electrons  is  2  cm  and  that  of  0+ 
ions  4  m  when  moving  perpendicular  to  a  0.  5  G  magnetic  field  with  0. 1  eV  energy. 
For  an  orbiting  cylinder  with  a  =  2.  5  m  and  ba  =  10  m  (the  approximate  size  of  a 
shuttle),  the  theory  would  hold  for  electrons  but  not  for  0+  ions.  So,  application 
is  made  to  a  cylinder  ten  or  a  hundred  times  larger.  (Earth  satellites  this  large 
are  anticipated  within  a  few  decades. ) 

Trajectories  of  the  guiding  center  ai  e  obtained  as  follows:  Eq.  (19)  and  E  * 
-70 are  used  in  Eq.  (22).  sin  0  is  set  equil  to  y/R,  and  the  result  is  equated  to 
the  limit  of  the  same  expression  »s  R  +  •  ,  Thl-  gives 


p(p  -  Cjy.)  = 


I  fp  ♦  +  - ! 

|  [  R2(b  +  R)2JR  R2  (1  + 

+  A.  zi\[P  — _ _ I2 

\  R2/[  R2(b  +  R) 


f2 _ 

■  b/R)1-u 


(bTm]y  * CV2 


*  C1  P  "  -9 
1  R2 


where  ym  is  the  limit  of  y  as  R  -  • ,  that  is,  the  impact  parameter; 


_  2a  _  2eBa 
C1  ~  Rg  ~  mVQ  ’ 


C2  caV0B’ 


Eq.  (23)  may  be  written  as 


Dy  +  Fy  +  G  =  0, 


& 

$ 

I* 

‘v* 


where  D,  F,  and  G  are  functions  of  R  and  are  independent  of  y.  The  solution  for 
y  and 


x  =  t  /r2  -  y2 


(27] 


give  the  trajectory  parametric  in  R.  The  "+"  sign  in  Eq.  (27)  has  physical  mean¬ 
ing  only  for  trajectories  that  do  not  hit  the  cylinder.  For  <t>c  not  too  negative  for 
electrons  and  not  too  large  for  positive  ions,  minimum  and  maximum  values,  y 
and  y  ^  ,  of  y „  exist  such  that  trajectories  starting  at  large  -x  for  which 

y„m  <  y„  <  y^  hit  the  cylinder.  Trajectories  with  y*.  <  y  pass  the  cylin¬ 
der  with  y  oc,<  -1  (at  x  *  0),  and  those  with  y«,  >  y^^  pass  the  cylinder  with 
y„  >  1.  Figure  4  shows  the  two  0+  trajectories  for  the  parameter  values  in 
Table  1,  a  =  25  m,  and  yM  =  yo>m  and  y^.  For  -x  starting  large  and  yoorn<  yoo 
^ y^M'  x  ^  ®  for  all  parts  of  the  trajectory.  Analytical  means  of  obtaining  ywn 

or  y  arc  not  feasible  for  these  tvpes  of  trajector ies.  They  may  be  obtained  to 
°°M  2 

any  desired  accuracy,  however,  as  follows:  From  Kqs.  (23)  to  (26),  set  F  - 

41X1  -  0.  solve  for  ,  and  obtain  its  minimum  value  as  R  is  varied.  This 
gives  yooM  .  The  value  of  yoom  is  obtained  by  solving  Kq.  (23)  for  ,  setting 
y  -  -R,  and  finding  the  maximum  value  ofyooas  R  is  varied.  Mathemati¬ 
cally.  bounded  ion  trajectories  exist  in  the  forbidden  regions,  the  shaded  ones  in 
Figures  4  through  6.  Since  no  ion  emission  is  assumed,  no  ions  should  exist 
there,  so  they  are  ignored  here.  The  sign  of  yF  -  4DG  in  the  solution  for  y  de¬ 
pends  on  whether  a  bounded  or  unbounded  trajectory  is  being  considered  or,  if  an 
unbounded,  whether  R  is  increasing  or  decreasing. 


Table  1.  Values  of  Parameters  Used  Not 
Specified  Otherwise 


parameter 

value 

parameter 

value 

c 

1.777 

B 

0.5G 

Vo 

7  km/  sec 

*c 

-0.5  V 

a 

25  or 

n 

1011  m'3 

250  m 

T 

2000°  K 

b 

4 

11 


Figure  5  shows  the  same  as  Figure  4  except  that  $c  ■  the  maximum 

value  for  any  ions  to  hit  the  cylinder.  For  this  case,  y  *  y  ... 

00  m  00  M 

Figure  6  shows  the  same  as  Figure  4  except  that  the  cylinder  is  ten  times 
larger.  The  limiting  trajectories  are  seen  to  be  different,  simpler  types.  For 


Y 


Figure  4.  Limiting  0+  Trajectories  for  B  *  0.5G,  VQ  ■ 
7  km/sec,  <pc  =  -0.5  V,  b  =  4,  and  a  =  25  m 


them,  yoom  and  yraM  may  be  obtained  from  the  maximum  and  minimum  values  of 
the  solution  of  Eq.  (23)  for  with  R  *  1  and  -la?  y.  ^  y,..  a=  1.  For  y  * 
y^ m  and  this  set  of  parameters  ($c  =  -0.5  V)  but  varying  a.  the  transition  between 
the  two  types  of  trajectories  in  Figures  4  and  6  occurs  at  a  *  46.  6  m.  This  is  de¬ 
termined  by  differentiating  Eq.  (23)  with  respect  to  R,  equating  it  and  dy/dR  to 
zero,  setting  R  =  1,  and  solving  (a  cubic  equation)  for  a.  The  maximum  ya  of  the 
forbidden  region  is  seen  to  be  0.5  or  0.8Rg  [see  Eq.  (24)]  from  the  cylinder  for 
a  -  25  or  250  m,  respectively. 

For  the  simple  type  of  trajectories  in  Figure  6,  the  difference  of  the  two  lim¬ 
iting  values  of  yM  is 


°M 


-  y0 


(b 


b  [" (b  >  l)Uc 

+  1)  jcjl  L  b 


e 

lei 


2b  +  3 
b  +  1 


(28) 


This  holds  when  the  bracket  is  not  positive  for  e  >  0  and  not  negative  for  e  <1  0, 

From  this  and  Eq.  (25),  as  expected,  it  is  seen  that  this  difference  increases  for 

electrons  and  decreases  for  ions  as  increases.  A  >,  is  the  value  of  A  for 

c  cM  c 

which  the  bracket  in  Eq.  (28)  vanishes,  y  -  y  is  zero  for  A  I A  „«  -»  1. 

'“M  J»m  rcycM 

For  the  previously  used  parameters,  A  .*  27.7  V  for  a  *  25  m  and  A  -  277  V 

CM  CM 
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for  a  *  250  m.  For  ions  and  0C  >  *CM*  the  value  of  y,,,  which  divides  the  trajec¬ 
tories  passing  the  cylinder  with  y  >  1  from  those  with  y  <  -1  is  determined  as 

described  above  for  obtaining  the  value  of  y  Table  2  shows  some  values  of 

°°  m 

yoom*  y“M*  and  their  difference,  along  with  those  of  ylm  and  y1M<  the  corre¬ 
sponding  limiting  values  of  y  where  the  particle  hits  the  cylinder.  The  results  for 
b  =  4  are  seen  to  be  close  to  those  for  b  =  “  .  As  b  -*■  °°  ,  the  term  containing  in 
Eq.  (28)  vanishes,  so,  although  the  trajectories  are  affected  by  0c>  y«M  '  y«m 
and  thus  the  current  per  unit  length  to  the  side  are  independent  of  it.  The  values 
of  y“M  '  y“m  for  °+  and  —  0  are  seen  generally  to  be  much  less  than  2,  the 
minimum  value  without  magnetic  field.  This  indicates  that  the  magnetic  field  de¬ 
creases  the  particle  currents  to  the  cylinder  and  the  y  extent  of  the  wake. 

For  simplicity,  the  current  to  the  cylinder  is  determined  assuming  only  cur¬ 
rent  from  ambient  electrons,  I  ,  and  0+  ions,  I+.  The  current  due  to  each  specie 
is  the  sum  of  Is,  the  current  to  the  sides  (curved  surface),  and  Ie,  that  to  the  ends. 
The  steady  state  potential  is  the  value  making  the  net  current  zero.  The  current 
to  the  sides  is  given  by  Eq.  (15).  For  each  specie, 

ls  =  neV  2ba"  (y  -  yV  (29) 

s  o  °°M  ,oom' 

This  gives  the  ram  current  without  magnetic  field  and  0c  =  0  ify„M  "  y«>m  =  2.  1^ 

vs  0c  for  the  parameter  values  in  Table  1  with  a  =  250  m  is  shown  in  Figure  7. 

Its  shape  is  parabolic  since,  from  Eq.  (25),  Eq.  (28)  and  thus  I*  are  quadratic 
functions  of  0  ,  at  least  for  the  simple  type  of  trajectories.  For  electrons,  b  =  4 

c  -  _4  _ 

and  a  =  250  m:  I  is  ~  10  of  I  and  thus  is  omitted.  Since  the  thermal  speed, 

Eq.  (32),  for  electrons  is  ten  times  V  ,  the  electron  current  may  be  calculated  ne¬ 
glecting  the  cylinder's  motion.  Since  the  electron  gyroradius  r  ~  2  cm  so  that 
.  g  g 

a  rg  ~  12»500>  according  to  Rubinstein  and  Laframboise  (consistent  with  the  top 

six  numbers  in  the  last  column  of  Table  2),  the  electron  current  to  the  sides  is 
negligible,  and  the  currents  to  the  ends  are 


T+  _  O  2 

le  ne+  VT+  27ra  • 


and,  for  <pc  <  0, 


=  ne.vT-  2na2  exp[-e_  0  /(kT_)], 


where  k  is  Boltzmann's  constant  and 


-  •  >  *  ■  »  »  *  , 


t  n  at  right  means  *  10 


0  50  100  150  200  250  300 

<$C  (Volt) 

Figure  7.  The  0+  Current  to  the  Side  of  the  Cylinder  as  a  Function 

of  <fi  for  a  =  250  m 
c 


1  +  I 


and,  for  steady  state. 


I  =  0,  (34) 

where  I*  is  given  by  Eqs.  (28)  and  (29),  For  a  =  25  or  250  m  and  the  parameter 
values  in  Table  1,  <f>c  =  -0.26  or  -0.61V,  respectively,  in  steady  state.  This  is 
determined  by  solving  Eqs.  (33)  and  (34)  for  the  0 in  Eq.  (31).  It,  and  thus  Ig, 
are  almost  independent  of  the  <p  in  in  Eq.  (28)  for  |  0c  |  <  <  I  #cM  I  • 

Now  the  time  response  resulting  from  capacitance  is  calculated.  Using 
Eqs.  (28)  through  (33), 

S1  =  1  ;  c  =  al  (~a2  +  a3*c)2  +  a4  -  a5  exp  (a6*c>*  (35) 

where  q  is  the  charge  on  and  C  the  capacitance  of  the  cylinder,  the  a.'s  >  0  and 

are  independent  of  <p  and  t.  The  value  of  C  is  obtained  from  Smythe's  slightly 

modified  empirical  formula  for  b  =  4  and  a  =  250  m.  The  result  is  61.8  nf. 

<b  ,  the  steady  state  value  of  0  at  t  =  <*>  ,  is  obtained  by  setting  Eq.  (35)  to  0  and 
C“>  c 

solving  for  the  <Pc  in  the  exponential  (as  near  the  end  of  the  previous  paragraph). 
Eq.  (35)  could  not  be  integrated.  However,  since  |0CI  <<  l0cMl, 
a3  |0C  I  <  <  a2  (indicating  that  I+  is  nearly  constant).  So,  by  replacing  0c  with 
a.^coo(with  very  little  error),  Eq.  (35)  could  be  integrated,  giving 


0„ 


_ h _ 

1  -  (1  -  h)  exp  (-  a^aght/C)  * 


(36) 


assuming  0c  -  0  at  t  =  0,  where 

h=  [aj  <-a2  +  a3  ^co>)2  +  a4]/a5  =  exp  {a6*cJ'  (37) 

Figure  8  shows  as  a  function  of  t  for  0  =  -0.  61V  and  the  parameter  val- 

C  Coo  C  CD 

ues  in  Table  1  with  a  -  250  m.  This  is  equivalent  to  the  potential  falling  exponen¬ 
tially  from  zero  to  0c  ^  with  a  time  constant  increasing  from  TQ  at  t  =  0  to  Tx  at 
t  =  “  ,  where 


the  charging  time  resulting  from  capacitance  is  small  compared  to  particle  transit 
times. 

1.  CONCLUSIONS 

Simplifying  assumptions  lead  to  a  relatively  simple  analytic  theory  for  the 
steady  electrical  state  of  a  body  moving  through  a  magnetoplasma.  The  first  two 
orders  of  two-dimensional  guiding  center  theory  in  the  object  frame  are  used.  It 
leads  to  the  potential  satisfying  Laplace's  equation  wherever  all  of  the  streamlines 
come  from  far  upstream.  It  is  applied  to  conditions  consistent  with  a  long,  25  m 
or  250  m  radius,  conducting,  right,  circular  cylinder  at  a  few  hundred  kilometers 
altitude  with  its  axis  parallel  to  the  magnetic  field  and  moving  perpendicular  to  it. 
Hot  electrons  but  cold  ions  are  assumed.  The  potential  is  obtained  from  an  empir¬ 
ical  formula  for  a  charged  cylinder  moving  across  a  uniform  magnetic  field  in  a 

e\ 

vacuum.  Choosing  f  =  1  -  b/[R  (b  -  1  +  R)]  would  make  Eqs.  (17)  and  (19)  slightly 
more  accurate.  No  charge  in  the  plasma  to  balance  that  on  the  cylinder  is  assumed. 
The  potential  field,  particle  trajectories,  and  currents  are  calculated.  The  mag¬ 
netic  field  decreases  the  ion  current  significantly.  Because  of  their  light  mass, 
nearly  all  electrons  that  approach  the  curved  surface  go  around  the  cylinder  with¬ 
out  hitting  it,  so  the  electron  current  to  the  sides  is  negligible.  The  currents  to 
the  ends  of  the  cylinder  are  the  same  as  without  any  magnetic  field. 

"Forbidden"  regions  next  to  the  cylinder  are  found  which  upstream  ions  do  not 
reach.  The  potential  caused  by  charge  density  in  these  regions,  ram  and  wake, 
and  its  effect  on  the  polarization  charge  are  neglected.  If  they  were  not  and  if 
charge  balance  were  assumed,  the  forbidden  region  at  y  >  -0.  5,  essentially  an 
electron  sheath,  would  be  smaller,  and  the  one  near  y  =  -1  would  probably  be  re¬ 
placed  by  a  positive  ion  sheath.  Since  these  neglected  effects  are  important,  the 
results  may  be  regarded  as  a  zeroth  order  approximation.  The  results  lead  to  tra¬ 
jectories  on  the  wake  side  being  symmetric  to  those  on  the  ram  side,  except  that 
the  wake  side  lacks  those  that  hit  the  cylinder.  (The  charge  in  this  region  of  the 
wake  is  also  neglected.  )  Thus,  there  is  no  particle  impact  on  the  wake  side  and 
no  ion  focusing. 

For  the  cylinder  charged  to  -0.5  V,  the  trajectories  are  simple  for  radii 
greater  than  47  m.  For  them,  the  ion  current  to  the  sides  increases  quadratically 
with  decreasing  cylinder  potential.  The  floating  potential  is  -0.  26  or  -0.  61  Y  for 
the  radius  equal  to  25  or  250  m,  respectively. 

The  charging  time  constant  resulting  from  the  capacitance  of  the  cylinder  is 
inversely  proportional  to  its  linear  dimensions.  For  a  cylinder  with  a  250  m  ra¬ 
dius,  simple  theory  results  in  less  than  a  nanosecond  time  constant.  This  is  neg- 
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ligible  compared  to  transit  times,  which  must  therefore  govern  the  charge -up 
time. 
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This  treatment  is  the  thick  sheath  approximation.  The  results  of  Smetana 
indicate  that,  at  least  for  no  magnetic  field,  less  than  10  percent  error  is  made  in 
the  attracted  species  flux  for  |£|  <  4kT/|e|.  The  thickness  of  the  forbidden  re¬ 
gion  at  y  >  -0.  5  is  the  order  of  R^  [  defined  in  Eq.  (24)].  Since  a  Laplacian  po¬ 
tential  is  used  throughout  the  whole  plasma,  the  effect  of  Debye  length  is  omitted. 

The  effect  of  charge  in  the  (much  thinner)  forbidden  regions  would  be  largely 
cancelled  by  balancing  the  charge  in  the  plasma  with  the  charge  on  the  cylinder. 

A  comparison  with  the  results  of  including  these  effects  in  two  or  more  iterations 
would  be  interesting.  This  would  indicate  the  accuracy  of  the  solution  here,  give 
a  better  solution,  and  indicate  how  rapidly  it  converges  to  the  latter.  The  larger 
cylinder  treated  here  has  a  radius  the  same  order  of  magnitude  as  the  mean  free 
path.  So,  either  collisions  should  be  taken  into  account  or  a  smaller  size  should 
be  chosen. 

Treatment  of  other  body  orientations,  shapes,  and  orbits  would  be  of  interest. 
Generally,  this  would  mean  attempting  to  extend  this  treatment  to  three  dimen¬ 
sions.  Since  the  velocity  of  a  satellite  at  200-  or  300 -km  altitude  is  only  an  order 
of  magnitude  or  less  larger  than  the  ion  thermal  motion,  including  the  latter  in 
the  treatment  would  be  desirable.  Waves  and  turbulence,  especially  for  large, 
irregular  bodies  such  as  the  shuttle,  may  exist  and  thus  also  should  be  considered. 
The  effect  of  the  finite  ion  (but  not  usually  electron)  gyroradius  must  be  included 
for  bodies  as  small  as  or  smaller  than  the  shuttle.  For  the  magnitude  of  cylinder 
charge  small,  the  electric  field  vanishes  along  two  lines  on  the  cylinder  parallel 
to  the  axis,  one  on  the  ram  and  one  on  the  wake  side.  For  the  magnitude  of  this 
charge  large,  the  electric  field  vanishes  on  a  line  away  from  the  cylinder  parallel 
to  the  axis  at  x  =  0.  Thus,  the  guiding  center  approximation  breaks  down  in  the 
vicinity  of  this  (these)  line(s).  Calculations  of  trajectories  that  come  close  to  this 
vicinity  may  thus  be  in  error.  For  comparison,  a  calculation  of  the  case  treated 
here  but  with  no  magnetic  field  would  be  interesting. 

A  calculation  using  the  POLAR  code  is  being  done  for  comparison.  Prelimi¬ 
nary  results  with  an  octagonal  cylinder  200  m  long  and  50  m  across  at  -0.26  V  po¬ 
tential  include  0.26  A  ion  and  0.03  A  net  current  compared  with  0.96  and  0  A, 
respectively,  here.  POLAR  gives  no  forbidden  region  on  the  ram  side  of  radial 
extent  larger  than  —  10  m. 


19.  Smetana,  F.O.  (1963)  On  the  current  collected  by  a  charged  circular  cylinder 
immersed  in  a  two-dimensional  rarefied  plasma  stream,  in  Rarefied  Gas 
Dynamics  (Vol.  II),  L.A.  Laurmann,  Ed.,  pp.  65-91. 
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